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The presence of a phase transition in a finite system can be deduced, together with its order, 
from the shape of the distribution of the order parameter. This issue has been extensively studied 
in multifragmentation experiments, with results that do not appear fully consistent. In this paper 
we discuss the effect of the statistical ensemble or sorting conditions on the shape of fragment 
distributions, and propose a new method, which can be easily implemented experimentally, to 
discriminate between different fragmentation scenarii. This method, based on a reweighting of the 
measured distribution to account for the experimental constraints linked to the energy deposit, is 
tested on different simple models, and appears to provide a powerful discrimination. 
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O ' I- INTRODUCTION 

(N 

?H I At the transition point of a first order phase transition, the distribution of the order parameter in the corresponding 
^ ^ finite system presents a characteristic bimodal behavior in the canonical or grancanonical ensemble [3, [3, 0|. The 
bimodaUty comes from an anomalous convexity of the underlying microcanonical entropy. It physically corresponds 
to the simultaneous presence of two different classes of events, which at the thermodynamic limit leads to phase 
' coexistence. Such behavior is very different from the signal of a second order phase transition, or more generally a 
I , critical behavior. In this latter case the order parameter distribution is generally wide but always monomodal, it 
J>. fulfills specific scaling properties with increasing system size, and presents non-gaussian tails 0. 

fSJ ] In the case of nuclear multi-fragmentation, the most natural observable to analyze as a potential order parameter 
' is the size of the heaviest cluster produced in each collision event. Indeed this observable is known to provide an 
■ order parameter for a large class of transitions or critical phenomena involving complex clusters, from percolation to 
gelation, from reversible to irreversible aggregation^. Moreover it is reasonable to believe that the largest cluster 
size is always well correlated to the total energy deposit; this means that it will take different values in the two 
, phases of any transition involving a finite latent heat, thus serving as an order parameter. Many different indications 
' exist (^,01 that multi- fragmentation may be the finite size and possibly out-of-equilibrium realization of the liquid-gas 
''pi I phase transition of infinite nuclear matter. Also in this case the heaviest cluster size is expected to play the role of 
' , ■ an order parameter, because of its good correlation with the system density, as shown by the numerical simulation of 
O I the liquid-gas transition[4|. 

^ These ideas have been recently applied by the INDRA collaboration to different intermediate-energy heavy-ion 

^ ' collisions data sets, with the aim of tracking the multi-fragmentation phase transition and possibly determining its 
J> , order [1, [l^, |Tl|, [l^, [l^, A systematic analysis of different data from central collisions reveals that the 

JJ"! ■ heaviest fragment distributions are never bimodal, neither they present the nonjgaussian tails that would allow to 
rS ' identify a critical phenomenon. The situation is different in peripheral collisions [9| . In this case, the quasi-projectile 
' heaviest residues are sorted in centrality bins measured from the transverse energy of light charged particles detected 
on the quasi-target side. The resulting distributions clearly show two different event families when plotted against 
the asymmetry between the two heaviest fragments Zasy = {Zmax — Z second) /{Z max + Z second) , and do not fulfill 
the scaling properties shown by central collisions. The widths of the two partial distributions are too large for a 
bimodality to appear in the projection over the Z„iax axis, but a clear bimodality is seen in the closely correlated 
Zasy distribution. 

The characteristics and order of the nuclear multi-fragmentation transition are thus not completely clear. 

It is certainly possible that two very different fragment formation mechanisms act in central and peripheral 
collisions [15]. It is however also important to note that the sorting conditions are not the same in the different 
collisions, and the shape of the distributions is obviously inffuenced by the constraints imposed by the sorting. 

In this paper we analyze the effect of the sorting, and propose a new method to detect a possible first order 
phase transition independent of the mechanism of energy deposit, and without the artificial assumption of a heat 
bath. Indeed the dependence on the sorting is a manifestation of the statistical mechanics concept of ensemble 
inequivalence, as we develop in detail below. Working out the relations between the different ensembles will allow us 
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to predict the transformation to apply to the bimodahty signal, in order to account for the experimental constraints. 



II. BIMODALITIES AND ENSEMBLES 




FIG. 1: Event distributions in the isobar lattice gas model. A system of 216 particles is considered at a subcritical pressure and 
a temperature close to the transition temperature. Central panel: canonical distribution in the largest cluster size and energy 
plane. The left and right figures give the projection of this correlation over the largest cluster size axis. Left: canonical ensemble; 
right: narrox constraint on the total energy as indicated by the two vertical lines of the central panel. Energy is given m units 
of the closest neighbors coupling. 



Observables distributions depend on the criterium used to collect the events, i.e. on the statistical ensemble. 
Concerning the issue of bimodality, a system experiencing a first order phase transition presents a bimodal distribution 
of its order parameter if and only if the associated observable is constrained only in average through a conjugated 
Lagrange parameter, and it is not fixed by a conservation lawfsl. [iGj. 

For example, the energy distribution of a finite systern, presenting in the thermodynamic limit a finite-latent-heat 
phase transition, is bimodal in the canonical ensemble 0, [l3]- The same is true for any other observable Z presenting 
a non-zero correlation with the total energy. If the correlation is loose and the system is small, each of the two peaks 
may have a large width, which may extend over values characteristic of the other phase, and mask the bimodality. 
However, the distance between the two peaks Zi and Z2 scales linearly with the number of particles: {Zi — Z2) oc N, 
while the width of each peak scales only with its square-root: az oc v^, meaning that a bimodal behavior will be 
recovered in the Z observable, for sufficiently large systems. 

In the microcanonical ensemble, the energy distribution is by definition a (5-function, and obviously cannot be 
bimodal. The distribution of any other observable Z will in general present a finite width, but if the correlation 
with the energy is strong, the bimodality will be lost. These general arguments are applied to the liquid-gas phase 
transition in Fig[Tl which shows the event distribution of the Lattice Gas Model in the canonical isobar ensembleQ 
at the transition temperature. This model presents a transition with finite latent heat between a high density liquid 
and a low density gas. The order parameter of the liquid-gas transition is known too be one-dimensional, meaning 
that energy and density, or equivalently volume 1^ oc , have to be correlated. 

This intuitive statement can be formalized in the framework of quantum statistical mechanics pisl [l9j. Indeed the 
result of the statistical measurement of any observable Z, associated to the operator Z is 



< Z >^= TtZD, (1) 

where D is the density matrix and Tr means the trace over the quantum Pock or Hilbert space. In the space of 
Hermitian matrices, the trace provides a scalar product 



<< Z\\D >>= TrZD. 



(2) 
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It is then possible to define an orthonormal basis of Hermitian operators {Zi} in the observables space, and to 
interpret the measurement < Zi >^ as the projection of the density matrix over the corresponding axis. Statistically 
independent observations can thus be associated to orthogonal observables in the Liouville space. 

As a practical consequence, in the liquid-gas phase transition the distribution at the transition temperature must be 
bimodal both in the volume and in the energy direction^. As anticipated in the introduction, the size of the largest 
cluster is also correlated to both volume and energy, and as such it can be used as an equivalent order parameter. 
This is shown in the left side of Fig[TJ at the transition temperature, the distribution of the largest-cluster size shows 
two peaks separated by a region of convexity. These observations are done in the canonical ensemble. If energy is 
constrained as in the microcanonical ensemble, the two phases are not accessible in the same ensemble of events any 
more, and bimodality is lost. This is not only true for the exact energy conservation implied by the microcanonical 
ensemble, but for any strong constraint on the energy. This is shown in the right side of Fig[TJ the distribution issued 
from a narrow energy constraint on the canonical event set is a normal distribution. 



III. BIMODALITIES AND CONVEXITY ANOMALIES 



This behavior is a consequence of ensemble inequivalence, which has been addressed by numerous works in the 
literature ^2(?|. In particular the presence of a conservation constraint on an order parameter has been reported to 
modify different phase transition observables, making a first order transition look like a continuous onep^. [2ll[2^. 

The prototype of an ensemble constraining the order parameter is the microcanonical ensemble. There, the phase 
transition can be univocally recognized studying the curvature properties of the density of states [23|. Indeed from 
the standard definition of the canonical ensemble 

pp{E)^W{E)eM-PE)Z-' (3) 

we can immediately see that a bimodality in the canonical energy distribution is exactly equivalent to a convex intruder 
in the microcanonical entropy S — log W^, which leads to the well known microcanonical negative heat capacity ^4, ,2^]. 
Let us now consider the case of a second observable Z. If both E and Z are order parameters, and the transition is 
first order, the two-dimensional probability distribution p{E, Z) should be bimodal in both the E and the Z direction 
within the ensemble where the observables are fixed by the conjugated Lagrange multipliers /3, A 0]: 

ppiE, Z) = W{E, Z) eM-PE - \Z)Z^\ (4) 

All conservation laws on other variables are implicitly accounted in the definition of the density of states W . For 
instance, if Z represents the largest cluster size and the total system size is Ztoti W reads 

W{E, Z) = Wz.jE, Z)^Y. - £^<"^)'5(^ - Z^^^^)5{Ztot ~ Z^ul) (5) 

(«) 

where the sum extends over the system microstates. The search for bimodalities can only be done in this (extended) 
canonical ensemble, and is exactly equivalent to the study of the curvature matrix of the entropy in the two-dimensional 
observable space 




If this curvature matrix has two positive eigen- values, this means that and Z and E are associated to two independent 
order parameters. In the more physical case of a one-dimensional order parameter, only one eigen-value is positive, 
and the associated eigen- vector can be taken as the '"best"' order parameter. It is the linear combination of the E 
and Z observables, which gives the best separation of the two phases in the two-dimensional spacefS]. 

In the physical case of nuclear multi-fragmentation experimentally studied through nuclear collisions, the distribu- 
tion of the deposited energy crucially depends on the entrance channel dynamics and data selection criteria. In the 
case of quasi-projectile events selected in heavy ion collisions, the energy distribution is very large, and is determined 
by the impact parameter geometry and dissipation dynamics. If events are sorted in centrality bins, the distribution 
is centered on a well defined value given by the average dissipation at the considered impact parameter, but the 
distribution has a finite width that depends in a non-controlled way on the selection criteria. The statistical ensemble 
describing multi-fragmentation data is thus neither canonical nor microcanonical, and should rather be described in 
terms of the gaussian ensemble [25|l, which gives a continuous interpolation between canonical and microcanonical. If 
E and Z have a non zero correlation as we expect , W{E^ Z) ^ We{E)Wz{Z), the distribution of energy will affect 
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also the distribution of and the concavity of the Z distribution will not be univocally linked to the concavity of 
the entropy. 

Concerning the explicit constraint A on there is no reason to believe that the collision dynamics or the data 
treatment induces a specific constraint on the size of the largest cluster, other than the total mass and charge 
conservation, which are already implemented in the definition of the state density. This means that we will consider 
in the following A = 0. 



IV. RELYING OBSERVABLES DISTRIBUTIONS TO THE UNDERLYING ENTROPY 



In the last section we have argued that, in the absence of a canonical sorting, there is no one-to-one correspondence 
between the shape of the probability distribution and the phase transition properties of the system implied by its 
density of states. This is not only true for the energy, but also for all other observables that present a non-zero 
correlation with the energy. Unfortunately, if the order parameter is one-dimensional, all observables that can be 
proposed to study bimodality (charge of the heaviest cluster, asymmetry, etc..) must be correlated with the energy. 
This means that the presence (or absence) of the bimodality signal may depend on the experimental sorting conditions. 

It is however important to note that if the energy distribution cannot be experimentally controlled, it can be 
- at least approximately - a-posteriori measured. This means that it is possible to unfold from the experimental 
distribution the contribution of the entropy, giving the phase properties of the system, and the contribution of the 
energy distribution, which depends on the collision dynamics. Indeed, as long as no explicit bias acts on the Z 
variable, the experimental distribution can be calculated from the canonical one eq.(|4]) by a simple reweighting of the 
probabilities associated to each deposited energy 

Pe.p(i?,^)-P,(^^>^)^^ (7) 

where PexpiE) is the measured energy distribution. 

In the following, we explicit the effect of this experimental bias in the two cases of interest: an entropy curvature 
matrix with two negative eigen-values, corresponding to the absence of a phase transition, and an entropy presenting 
one direction of convexity (first order phase transition) . 



A. The monomodal case 




FIG. 2: Canonical distributions in the plane of the total thermal energy, normalized to the total available energy, and the size 
of the largest cluster, normalized to the total system size, tn the (double) saddle point approximation. Right part: case of a first 
order phase transition eg. f20\) . Left part: without transition, eg. {Pp. The corresponding parameters are given in the text. 

Let us consider again the density of states W{E, Z) — ex-p{S{E, Z)) in the two-dimensional observable space defined 
by the total energy E and the largest cluster size Z . In the absence of a phase transition, this function is concave 
everywhere. To evaluate the partition sum at a temperature I3~^ we can make a saddle point approximation 

S{E, Z) = Sq + I3{E - Bp) - cii{E - Epf - C22{Z - Zpf + cu{E - Ef,){Z ~ Zp) (8) 
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where Ep^Zp are the most probable values of E and Z at a temperature /3 ^, and C = {cij} is the entropy curvature 
matrix eq.®. In this approximation the event distribution is a gaussian 



PpiE,Z) = 



1 



1 



: exp 



1 



where x = {E — Ep, Z — Zp)^ and the variance- covariance matrix E is related to the curvature matrix by 

1 d'^S. 

1 d^S , 
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(1 - p^)(Je(^z 



Cl2 = C21 



2dEdZ 



(9) 

(10) 

(11) 
(12) 



An example is given in the left part of Figure [H To fix the ideas we have chosen Ep = 0.35£^(ot: Zp = QAbZtot, 
ge — 0-^Etot, <^z — 0-^5Ztot, P = —0.6, where Etot and Ztot are the total available energy and size, respectively. The 
inclusive probabilities associated to this distribution, or marginal distributions, are simple gaussians 



1 



V2w 
1 



exp 



v/2^ 



exp 



{E ~ Epf 

{Z-Zpf 
2a| 



(13) 

(14) 
(15) 



Deviations from a perfect gaussian are in principle possible because of the effect of the boundary, which limits the 
possible values of E and Z due to the conservation laws. These effects are however very small in all physical cases of 
interest, and the marginal distributions are largely independent of the correlation coefhcient. 




FIG. 3: Marginal distribution for the monomodal case represented in the left part of Figure\^ Left: gaussian ensemble eq.^. 
Right: reweighted distribution eg. \24^ . The correlation coefficient is varied from p = —0.15 (black lines) to p — —0.85 (lightest 
grey). The other parameters are taken as for Figure\^ 

As we have stressed in the previous section, in the experimental case, because of the absence of an external heat 
bath, there is no reason to believe that the event distribution can be described by eq. ([§])• Even if the partitions may 
fully explore the available phase space, the dynamics of the entrance channel plus the detection and sorting constraints 
will lead to an energy distribution which is not uniquely characterized by its average value < E >. When data are 
sorted in centrality bins, the energy distributions as measured from calorimetric algorithms are generally close to 
gaussians. This can be easily understood in the framework of information theorvp^. where such a situation can be 
dealt with introducing the energy variance as a second constraint, which naturally leads to the gaussian ensemble [2^. 
Therefore we will take as a typical shape for pexp(£'j Z) the form ^ where Pexp{E) is a gaussian of average i?exp and 
variance Ceajp- This choice is only done to present some realistic numerical applications, and is readily extended to 
more general situations. 
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Eq.® and eq.([7|) have the same functional dependence, with an important difference though. In the case of eq.® 
the first two moments of the energy distribution eq. (|13p are directly related to the canonical partition sum In by 



dliiZ. 



13 



dp 
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and, in the saddle point approximation, to the microcanonical temperature T and heat capacity 



T- 



dS{E,Z) 
dE ' 



d^S{E,Z) 
dE^ ' 



(16) 
(17) 

(18) 
(19) 



On the other hand, the first two moments of the energy distribution ([7]) are external constraints imposed by the 
collision dynamics, that have a priori no direct connection with thermodynamics. Since the folding of gaussians is still 
a gaussian, the corresponding Z distributions, shown for the model case discussed above in the left part of Figure [3l 
are also of gaussian shape. The width of the distribution depends on the value of the correlation coefficient p (a lower 
value of the correlation leading to a larger distribution), and all distributions are normal. For the specific application 
of Figure [31 we have fixed E^xp = Ep, Uexp = o'b/2. 

These observations mean that, if the system does not present a phase transition, the non-canonical experimental 
sorting does not qualitatively bias the shape of the measured Z distribution. 



B. The bimodal case 




FIG. 4: Marginal distribution for the bimodal case represented in the right part of Figure\^ Left: gaussian ensemble eg. S22\) . 
Right: reweighted distribution eq. \24^ . The correlation coefficient is varied from p = —0.15 (black lines) to p — —0.85 (lightest 
grey). The other parameters are taken as for Figure\^ 

In the presence of a first order phase transition with a finite latent heat, the saddle point approximation eq.® 
fails [23t. If the system is not too small, a double saddle point approximation around the two maxima {Ei, Zi), {Eg, Zg) 
can however be employed leading to a double humped distribution 

ME,Z) . A^.^exp +A^.^exp [-'-.-g^g^) (20) 

where xl = (E — Ei,Z — Zi),i — l,g, (Sg) represent the variance-covariance matrix evaluated at the liquid (gas) 
solution, and Ni{l3),Ng{f3) are the proportions of the two phases , with Ni defSg = Ng^/detY^i at the transition 
temperature f3t . An example of a distribution at j3t is given in the right part of Figure [21 To fix the ideas we have 
chosen El = 0.2Etou Zi = O.TZt^u <Jei = O.UEtot, <Jzi = 0.15Zt,uEg = O.bEt^u Zg - O.lSZ^^t, asg = O-ObEtou 

azg = Q.lZtot, Pi = Pg- 
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If the distance in the Z direction between the two phases Zi — Zg is sufficiently large respect to the associated 
variances ozi^ '^Zg^ as it is the case in Fig|5J the presence of the transition can then be inferred looking for the 
temperature interval associated to a bimodal shape of the Z distribution 

pp{Z) = AT, J£^exp ( ) + A^,,/:^exp (21) 

If we do not dispose of a canonical sample however the situation is not so simple. Let us consider the case where the 
experimental energy distribution belongs to the gaussian ensemble with an average E^^j, and a variance (J^^^p which 
can be experimentally measured for each given sorting. Then the measured Z distribution reads (taking pi = pg): 

Pe.,{Z)^ dEpp,{E,Z) ) (22) 

JEr^^r. expf- ^'l-f'^ ) + exp(- <^'fi^ ) 

which may look close to gaussian even in the phase transition region. In particular if the energy distribution is not 
very large, the sample corresponding to the transition energy Et (i.e. the energy corresponding to the minimum 
probability of eq. ([20|l ) will not explore the energy domains corresponding to the liquid E k, Ei ± aEi and to the gas 
phase E K Eg ±aEg- If Z is well correlated to E as we expect, then the only size partitions explored at the transition 
energy will be the intermediate ones between the liquid and the gas, leading to a normal Z distribution. This can be 
clearly seen in the left part of FiglH Only for very low values of the correlation coefficient the presence of the two 
phases can be recognized from the Z distribution, where the minimal value for p depends on the width of the inclusive 
energy distribution in the specific experimental situation under study. 

V. REWEIGHTING THE ENERGY DISTRIBUTION 

In order to recover the information on the concavity of the entropy from the measured distribution eq. (j22[) , we have 
to get rid of the dominant energy dependence PexpiE). To this aim we can reweight the total {E, Z) distribution in 
such a way that the total energy distribution is a constant between Emm and Emax, with Emin < E < E^ax, and 
Emim Emax are chosen large enough, such that the spanned energy interval contains the two phases. This procedure 
corresponds to introducing a weight for each E bin 

w{E) = ^ dZ p,,piE, Z)^ (23) 

and generates a new weighted distribution which is a direct measure of the density of states: 

pUE,Z)^'-^ = ^^. (24) 

In the limiting case of a perfect (anti)correlation between E and Z, the concavity anomaly in the two-dimensional 
entropy S(E, Z) is solely due to the convex intruder in the function S{E). In this case the order parameter is aligned 
with the E direction, and it is clear that the reweighted distribution eq. (|24p will be a normal distribution in spite of 
the presence of a first order phase transition. However, in the case of the liquid gas phase transition we have good 
reasons to expect [26j the order parameter in the enthalpy H — E +pV direction. Since Z is correlated to the volume 
as well as to the energy, then we may expect a non zero component of the order parameter along the Z axis, and an 
imperfect correlation between Z and E. In this case an anomaly in S{E, Z) may persist even after subtraction of the 
entropy S{E), and this anomaly will be shown by a residual bimodality in the distribution (j24p . 

A. Application to model cases 

The projections of p^ in the Z direction are shown in a perfectly monomodal and perfectly bimodal case, in the 
right parts of FigslS] and |4] respectively. If the underlying entropy presents a convex intruder (FigH]), the bimodality 
persits in the reweighted distribution for all considered values of the correlation coefficient. The case of a normal 
distribution is explored in FigO the distribution is still approximately gaussian if the correlation coefficient is small. 
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FIG. 5: Distributions of the size of the largest cluster for a system of N=216 particles within the isobar lattice gas model at 
subcritical pressure. Left: canonical calculations at temperatures T = 3.72,3.74,3.76,3.78,3.80 MeV close to the transition 
temperature (thin lines), and their sum (thick line) obtained weighting evenly the different temperatures. Right: microcanonical 
distributions at energies spanning the coexistence region E/N = 0.75, 1.25, 1.75, 2.25, 2.75 MeV (thin lines), and their sum (thick 
line) obtained weighting evenly the different energies. The closest neighbors coupling is fixed to e = 5.5 MeV. 



while if p is close to —1 it tends to a constant, but in no case the reweighting procedure applied to the monomodal 
distribution creates a spurious bimodality. 

Another example of the effect of the reweighting procedure is given in Fig[51 The right part of Fig[5] shows 
microcanonical calculations within the lattice gas modelQ]. Several distributions of the largest cluster size obtained 
for different total energies are displayed. In agreement with Fig[Tl these distributions are never bimodal, even if the 
total energy interval is chosen such that the first order phase transition is crossed by the calculation. The bimodality 
is recovered if these distributions are summed up with equal weight following eq. ([M|) (thick line in FiglS]). The detailed 
shape of the reweighted distribution obviously depends on the energy interval used, a wider interval leading to an 
increased distance in energy between the maximum and the minimum. 

For comparison, canonical distributions at different temperatures close to the transition temperature are displayed 
in the left part of the same figure. Because of the finite system size, the distributions are bimodal over a narrow but 
finite temperature range, while the transition temperature can be identified from the requirement that the two peaks 
have approximately the same height 3. We can see that the canonical distribution close to the transition point is 
wider than the reweighted microcanonical one, and its concavity anomaly is more pronounced, however the shapes of 
the two distributions are very close. The width of the transition interval (i.e. the distance between the two peaks) as 
estimated from microcanonical calculations is a lower limit of the physical transition interval, and the quality of this 
estimation depends on the energy width of the available sample. 

In the experimental case, due to the incomplete detection and the imperfect emission source selection, only an 
approximation of the deposited energy can be measured via measurable observables like transverse energy or calori- 
metric energy. We have checked that allowing a finite energy width does not change the result: an equal reweighting 
of gaussian ensembles leads to a distribution virtually identical to the one displayed in the right panel of FigO An 
equal reweighting of canonical distributions is also shown in Fig|5l the resulting distribution is almost superposable 
to the distribution at the transition temperature. 

These examples suggest that, if the microcanonical entropy presents a convexity, it may be possible to recognize 
it from the experimental data through a reweighting procedure, provided we dispose of an experimental sample 
sufficiently large in deposited energy. 

A more challenging situation is given by the presence of a critical phenomenon. In this case a variation of the control 
parameter leads again the system to pass from one phase to another, but this transition is continuous and does not 
corres pon d to a convex microcanonical entropy. Fig[5] shows calculations performed within the bond-percolation 
model|27|. This model presents a critical behavior at a value Pc of the bond-breaking probability. As the Lattice 
gas model can be considered as a prototype for phase coexistence and first-order phase transitions, similarly the 
percolation model is the simplest realization of a second order phase transition. The percolation distribution close to 
the critical point is very close to the microcanonical lattice gas distribution in the middle of the coexistence region 
of the first-order transition. The consequence of this well-known featu re [21I [2^ is that it is still not clear whether 
nuclear-multifragmentation is better described as a critical phenomenon [28!| or a first-order phase transition(29l. [soj . 

Figureinishows that also in this case the presence of a phase transition is detected by the reweighting procedure: a flat 




FIG. 6: Calculations in the bond three-dimensional percolation model for a system of 216 particles. Thin lines: distribution of 
the size of the largest cluster for equally spaced values of the bond-breaking probability from p — .16 to p = .56 including the 
critical value pc- Thick line: summed distribution obtained weighting evenly the different p's. 



distribution of bond-breaking probabilities leads to a global distribution qualitatively very similar to the distribution 
at the critical point if the p's interval is narrow, but if we have a sample including p's sufficiently far away from the 
critical probability as in Figure [SI two peaks appear corresponding to the ordered {p « 0) and desordered (p « 1) 
phase (thick line). 

From Figure [5] and Figure [H] we can see that in both models the reweighted distribution show the presence of the 
two phases, with an important difference though. If the lattice-gas calculation shows a clear convexity between the 
two phases, associated to the convex intruder in the microcanonical entropy, this is not the case for the percolation 
distribution which, in the critical mass region, largely keeps the convexity of the critical point, i.e. it is essentially 
flat. This statement is quantified in Figure [71 which shows the numerically calculated first and second derivative of 
the reweighted distributions in the two models. Only the model associated to a first order phase transition exhibits 
a region of positive convexity. 

This result can be easily understood. Indeed from eq. we can identify the (inverse) curvature of the reweighted 
distribution with a (generalized) susceptibility: 



(25) 



The back-bending shown by the lattice gas calculation in Fig[7| has then the same physical meaning as the well- 
known back-bending of the microcanonical caloric curve [2^, and allows to unambiguously identify a discontinuous 
(first order) phase transition. 



B. Application to experimental data 

In the previous section we have seen that, for different model cases that can be relevant for nuclear multi- 
fragmentation, the concavity properties of the reweighted distribution reflect the concavity of the underlying entropy. 

However we may wonder whether the different behaviors shown in Figs l3l4l5l6l will be distinguishable in realistic 
experimental cases, where distributions are often affected by important statistical as well as systematic error bars. 

Even more important, we may also ask whether the behaviors of the schematic models we have shown can be taken 
as general representatives of the presence or absence of a first order phase transition. From a mathematical point 
of view, the two-dimensional function logPu)(a;i, 2:^2) is concave everywhere if and only if all the eigenvalues of the 
curvature matrix c^^- = (S'(a;i,a;2) — S(xi)) /dxidxj, are negative. A necessary and sufficient condition that an 
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FIG. 7: Upper part: reweighted distributions of the size of the largest cluster for a system of N=216 particles with the lattice 
gas and the percolation model. Lower part: first (left) and second (right) derivatives of the distributions. 



hermitian matrix has negative eigenvalues is that the minors of the determinant are negative, i.e. 

detC < 
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which in our case becomes 



d'^S{E, Z) d^S{E) 



dE^ 



< 



<0 ; detC 



d^S{E) d'^S{E, Z) 



dE^ 



dZ^ 



< 



(26) 



for every {E, Z) value. In principle these conditions can be violated even if S{E, Z) and S{E) are concave everywhere. 
This means that we cannot a priori exclude a concave entropy even if we observe a bimodality in ■ 

In addition to that, only if a distribution is approximately symmetric we can visually judge, or numerically calculate, 
its convexity properties. In the experimental case, we may not dispose of an energy sample wide enough to construct 
a symmetric distribution around the transition value. It is clear for exemple that if in FigsO [6] we would sum up 
only low energy (respectively, low ps) samples, in both models a shoulder at low Z would be observed, and we would 
not be able to clearly discriminate between the first order and the second order scenario. This is probably the case 
in quasi-projectiles events, where the dynamics of the collision strongly favors low excitation energy deposit, and the 
desordered phase may barely, if ever, be attained [o']. 

To be conclusive about the presence of a first order phase transition we have therefore to show that the observed 
distribution is not compatible with a concave entropy, or in other words that a single saddle point approximation 
cannot explain the data as we have developed in section 4. 

In realistic experimental cases this can be tested through a test. Indeed in the absence of a phase transition 
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the underlying entropy is everywhere concave, a single saddle point approximation is reasonable, and the reweighted 
two-dimensional probability should be well described by: 



W{E,Z) _ 1 



exp 



'2(l-p2) [^J^ 



En 



Zn) 



(27) 



The parameters Eq, ZQ,aET<^ZT P should be fitted to the experimental reweighted distribution. A test should 
allow to exclude the adequacy of the fit. 

If fragmentation events are compatible with a critical phenomenon, both the single and the double saddle point 
approximation should fail. This means that in this scenario we should not be able to describe the data with the 
reweighted single gaussian eq . ((27|) . not with the reweighted double gaussian. An application of this method to 
INDRA data is in progress 



VI. CONCLUSIONS 

In this paper we have critically analyzed the bimodality observable proposed in J3] to track a first order phase 
transition, which has been extensively used to analyze multi-fragmentation data[l, [l^E, [lHj We have shown 
that the out-of-equilibrium dynamics of the entrance channel, as well as the different data sorting conditions, can bias 
the signal in an uncontrolled way, tending to suppress the bimodal behavior. Based on the correlation between the 
bimodality observable and the deposited energy, we propose a reweighting method on the measured distributions, with 
the aim of disentangling the dynamic dissipation properties from the entropy characteristics. This method applied 
on different schematic models, is seen to provide an excellent discrimination between the case of a first-order phase 
transition, a critical phenomenon and the absence of transition. A chi-square test is proposed for the experimental 
case to demonstrate/refute the existence of a convexity in the microcanonical entropy. 
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